home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_syr.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  3.5 KB  |  172 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_syr (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int uplo = 121;
  14.    int N = 1;
  15.    int lda = 1;
  16.    float alpha = 0.1f;
  17.    float A[] = { -0.291f };
  18.    float X[] = { 0.845f };
  19.    int incX = -1;
  20.    float A_expected[] = { -0.219597f };
  21.    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  22.    {
  23.      int i;
  24.      for (i = 0; i < 1; i++) {
  25.        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1402)");
  26.      }
  27.    };
  28.   };
  29.  
  30.  
  31.   {
  32.    int order = 101;
  33.    int uplo = 122;
  34.    int N = 1;
  35.    int lda = 1;
  36.    float alpha = 0.1f;
  37.    float A[] = { -0.291f };
  38.    float X[] = { 0.845f };
  39.    int incX = -1;
  40.    float A_expected[] = { -0.219597f };
  41.    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  42.    {
  43.      int i;
  44.      for (i = 0; i < 1; i++) {
  45.        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1403)");
  46.      }
  47.    };
  48.   };
  49.  
  50.  
  51.   {
  52.    int order = 102;
  53.    int uplo = 121;
  54.    int N = 1;
  55.    int lda = 1;
  56.    float alpha = 0.1f;
  57.    float A[] = { -0.291f };
  58.    float X[] = { 0.845f };
  59.    int incX = -1;
  60.    float A_expected[] = { -0.219597f };
  61.    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  62.    {
  63.      int i;
  64.      for (i = 0; i < 1; i++) {
  65.        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1404)");
  66.      }
  67.    };
  68.   };
  69.  
  70.  
  71.   {
  72.    int order = 102;
  73.    int uplo = 122;
  74.    int N = 1;
  75.    int lda = 1;
  76.    float alpha = 0.1f;
  77.    float A[] = { -0.291f };
  78.    float X[] = { 0.845f };
  79.    int incX = -1;
  80.    float A_expected[] = { -0.219597f };
  81.    cblas_ssyr(order, uplo, N, alpha, X, incX, A, lda);
  82.    {
  83.      int i;
  84.      for (i = 0; i < 1; i++) {
  85.        gsl_test_rel(A[i], A_expected[i], flteps, "ssyr(case 1405)");
  86.      }
  87.    };
  88.   };
  89.  
  90.  
  91.   {
  92.    int order = 101;
  93.    int uplo = 121;
  94.    int N = 1;
  95.    int lda = 1;
  96.    double alpha = -0.3;
  97.    double A[] = { -0.65 };
  98.    double X[] = { -0.891 };
  99.    int incX = -1;
  100.    double A_expected[] = { -0.8881643 };
  101.    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
  102.    {
  103.      int i;
  104.      for (i = 0; i < 1; i++) {
  105.        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1406)");
  106.      }
  107.    };
  108.   };
  109.  
  110.  
  111.   {
  112.    int order = 101;
  113.    int uplo = 122;
  114.    int N = 1;
  115.    int lda = 1;
  116.    double alpha = -0.3;
  117.    double A[] = { -0.65 };
  118.    double X[] = { -0.891 };
  119.    int incX = -1;
  120.    double A_expected[] = { -0.8881643 };
  121.    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
  122.    {
  123.      int i;
  124.      for (i = 0; i < 1; i++) {
  125.        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1407)");
  126.      }
  127.    };
  128.   };
  129.  
  130.  
  131.   {
  132.    int order = 102;
  133.    int uplo = 121;
  134.    int N = 1;
  135.    int lda = 1;
  136.    double alpha = -0.3;
  137.    double A[] = { -0.65 };
  138.    double X[] = { -0.891 };
  139.    int incX = -1;
  140.    double A_expected[] = { -0.8881643 };
  141.    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
  142.    {
  143.      int i;
  144.      for (i = 0; i < 1; i++) {
  145.        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1408)");
  146.      }
  147.    };
  148.   };
  149.  
  150.  
  151.   {
  152.    int order = 102;
  153.    int uplo = 122;
  154.    int N = 1;
  155.    int lda = 1;
  156.    double alpha = -0.3;
  157.    double A[] = { -0.65 };
  158.    double X[] = { -0.891 };
  159.    int incX = -1;
  160.    double A_expected[] = { -0.8881643 };
  161.    cblas_dsyr(order, uplo, N, alpha, X, incX, A, lda);
  162.    {
  163.      int i;
  164.      for (i = 0; i < 1; i++) {
  165.        gsl_test_rel(A[i], A_expected[i], dbleps, "dsyr(case 1409)");
  166.      }
  167.    };
  168.   };
  169.  
  170.  
  171. }
  172.